Precision polarimetry with real-time mitigation of optical-window birefringence 
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Optical-window birefringence is frequently a major obstacle in experiments measuring changes in 
the polarization state of light traversing a sample under investigation. It can contribute a signal 
indistinguishable from that due to the sample and complicate the analysis. Here, we explore a 
method to measure and compensate for the birefringence of an optical window using the reflection 
from the last optical surface before the sample. We demonstrate that this arrangement can cancel 
out false signals due to the optical-window birefringence-induced ellipticity drift to about 1%, for 
the values of total ellipticity less than 0.25 rad. 
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I. INTRODUCTION 

Kerr effect, measured in liquid helium [l[ and liquid 
nitrogen 0], is a useful method for non-contact mea- 
surement and monitoring of the electric fields in low- 
temperature experiments. In particular, we have pro- 
posed to use this technique to monitor the electric field 
set-up in the new search for the neutron electric dipole 
moment 0j. However, stress-induced birefringence of 
optical windows puts a limit on the sensitivity of this 
method. While constant stress due to the window mount- 
ing can be taken into account in modeling the signal 
by prior measurement of the optical-window birefrin- 
gence [4, 5] , or optically compensated in the experimental 
setup by use of retarder compensators [fHH, for exper- 
iments requiring greater sensitivity, the time-dependent 
stress due to temperature drifts cannot be treated us- 
ing the same methods, because the amount of retarda- 
tion due to the optical window at a given moment is not 
known. 

A careful experimental design, such as using low-stress 
window mounts and controlling the window tempera- 
ture, can minimize this birefringence-induced polariza- 
tion drift to as little as 3.5 x 10 -5 rad over the measure- 
ment time of 12 h [{|. But this is still not good enough for 
more sensitive experiments, such as measurement of Kerr 
effect in liquid helium. In a practical experimental setup, 
we expect Kerr effect signals of the order of 10~ 5 — 10~ 6 
rad LUl. In our laboratory, where low-tcmpcrature opti- 
cal access is implemented with a set of 1-in. -diameter, 
2-mm-thick, epoxy-mounted fused-silica windows, test- 
ing with 632 nm light showed that window birefringence 
introduces offsets in polarization rotation and ellipticity, 
the latter being defined as inverse tangent of ratio of the 
axes of the polarization ellipse, on the order of 0.1 rad, 
which also depend on the window temperature. This off- 
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set drifts within a range of about 10~ 4 rad over periods 
of hundreds of seconds, and this noise overwhelms the 
typical signal from the Kerr effect of liquid helium. 

In past experiments, with the exception of those where 
the residual birefringence in instruments is small enough 
to be ignored [lOj, various techniques have been used 
to distinguish the effects of the sample from those of 
the optical-window birefringence. In experiments with 
an isotropic sample, the incoming polarization of light 
can be rotated to distinguish the anisotropic feature of 
optical- window birefringence [ill ]. In other experiments, 
notably Kerr-effect experiments, modulation of the ap- 
plied electric field is used to distinguish the signal from 
background noise [l|. Neither of these methods can be 
used in the application of Kerr effect for monitoring elec- 
tric fields. Thus, we explore the feasibility of measuring 
and compensating for the optical-window birefringence 
drifts using the reflected light from the last optical sur- 
face before the sample. 



II. EXPERIMENTAL SETUP 

A typical experimental situation is simulated with two 
fused-silica parallelepipeds squeezed along the (approx- 
imately vertical) axis, making an approximate angle of 
7r/4 with the azimuth of the incident light polarization. 
These are labeled as "optical window" and "sample" in 
the schematic diagram of our experimental setup shown 
in Fig. Q] In the absence of stress on the optical win- 
dow, the maximum ellipticity that can be introduced by 
the sample on a single pass (depending on orientation) 
is equal to half of the birefringence-induced phase dif- 
ference, given by, (f> = 2nd(n s — n/)/A, where d is the 
thickness of the sample, A is the wavelength of the probe 
beam, and n s and rif are the indices of refraction along 
the slow and fast axes, respectively. We place a mirror 
just after the sample, so that the light double-passes the 
sample, and the entrance window is also the exit win- 
dow. This arrangement works for measuring the Kerr 
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effect since the effect of a small birefringence is additive 
for the double pass. The sample birefringence is what 
we want to measure, but this signal is obscured by the 
birefringence of the window. To mitigate this problem, 
we measure the window birefringence using the reflection 
from the back surface of the window. 

The laser light at 632 nm with intensity of 1 mW is 
produced by a Hitachi laser diode HL6320G mounted on 
a Thorlabs TCLDM9 mount, with laser current and tem- 
perature controlled by Thorlabs LDC 201ULN and TED 
200 controllers, respectively. The light is collimated to 
a beam width of about 1 mm diameter and is linearly 
polarized at 7r/4 rad from the optical axis of the sam- 
ple by passing through a Glan-Thompson polarizer. The 
beam is reflected back by a near-normal-incidence di- 
electric mirror placed after the sample to pass through 
the optical-window and the sample twice. Then it goes 
through a fused-silica photoelastic modulator (PEM) op- 
erating at 50 kHz, whose optical axis is aligned vertically, 
and a Wollaston-prism analyzer, which makes an angle of 
7r/4 rad with the optical axis of the PEM and, in this par- 
ticular setup, is crossed with the Glan-Thompson polar- 
izer. Additional dielectric mirrors were placed to ensure 
that all angles of reflection before the analyzer are less 
than 5° and to overcome space limitation for the mod- 
ulation polarimeter. For angles of incidence up to 20°, 
a commercial polarimeter, Thorlabs PAX5710VIS, was 
used to check that the dielectric mirrors used here do 
not alter the polarization state within the sensitivity of 
1 mrad for cllipticity and 0.1 mrad for polarization rota- 
tion. For lock-in detection, and for calibration purposes, 
we introduce light modulation at 700 Hz with a chopper 
wheel. The measurement and the reference beams are 
calibrated separately (as explained below) using the syn- 
chronously detected signals from each photodiode. All 
data are collected with a Lab VIEW program through 
the General Purpose Interface Bus (GPIB). 

The main experimental limitation is posed by the fact 
that the window birefringence is nonuniform. This can 
impart different ellipticity to the measurement beam and 
the reference beam (Fig. [TJ. To minimize the effect of 
this nonuniformity, we superimpose the two beams as 
they pass through the window. In order to be able to 
separate the measurement and reference beams, as well 
as the reflection from the front surface of the window, a 
wedged window is used (a wedge angle of 1° is sufficient 
in our setup). In our setup, we were able to reduce the 
distance between the two beams as they pass through 
the window to roughly one-tenth of the ~ 1 mm beam 
diameter. 

The photodiode signal with the modulation polarime- 
ter setup is given by, 
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I{t) = ^ {1 — cos(2e) sin(2a) cos (^4 sin w£) 
+ sin(2e) sin( A sin cot) } , 



(1) 



where e and a are the polarization ellipticity and azimuth 
of the light incident on the PEM, Iq is the intensity after 
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FIG. 1: Schematic diagram for the setup 



the polarizer, before the window and the sample, and 
K represents the reflection and absorption losses by the 
optical components before the photodiode. With a ~ 
7r/4 if no additional rotation is introduced by the sample 
and the optical window, A is the peak differential phase 
retardation introduced by the PEM, and uj ~ 2tt x 50 kHz 
is the PEM modulation frequency. The sine and cosine 
terms can be expanded in terms of Bessel functions of 
the first kind, multiplied by odd and even harmonics of 
w, respectively. In the present experiment, we detect the 
dc signal and the first harmonic signal. These two si gnal s 
are used for sensitive detection and calibration of e [12j • 



J(t) 
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{1 - cos(2e) sin(2a)J (A) 



+ 2sin(2e)Ji(A) sin(wi) H }. 



(2) 



Because the first-harmonic term depends only on e and 
is completely insensitive to a, this is a good polarimeter 
setup for measuring the light cllipticity. Furthermore, 
because 2.405 rad is a zero of the zeroth-order Bessel 
function of the first kind (Jo(2-405 rad) « 0), when the 
optical axis of the PEM forms the angle of 7r/4 rad with 
the analyzer, setting A = 2.405 rad ensures that the dc 
offset is independent of incoming light-polarization [l2| . 
This fact can be used to check and calibrate the PEM 
modulation amplitude to 2.405 rad by, for example, plac- 
ing a wave plate before the PEM and rotating it. PEM 
modulation amplitude is varied until the dc offset does 
not change with the rotation of the wave plate. In our 
setup, this method was used to calibrate the PEM re- 
tardation with approximately 1% accuracy, and the vari- 
ation of PEM amplitude within the separation between 
the two beams as they traverse the PEM is within this 
accuracy. For this special value of A — 2.405 rad, Eq. 
© simplifies to 



I(t) 
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{1 + 2J 1 (A) sin(2e) sin(wt)}. 



(3) 



The modulation polarimeter is thus calibrated for mea- 
suring ellipticity by dividing the intensity detected at the 
first harmonic of the PEM operating frequency, u>, by the 
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dc offset. Thus, e is given by 
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where A = 2.405 rad, and I\ and Idc are the intensity 
at the first harmonic of PEM operating frequency, mea- 
sured with a Signal Recovery 7265 lock-in amplifier, and 
the DC offset, measured at 700 Hz with a second Signal 
Recovery 7265 lock-in amplifier referenced to the chopper 
wheel, respectively. Before taking the ratio of two signals 
thus measured, we multiply the 50 kHz signal by an ad- 
ditional factor of 2 to account for reduction in intensity 
by 700 Hz chopping. Because the intensity signal itself 
contains the calibration information, this calibration can 
be done for both beams separately and, with additional 
lock-in amplifiers, simultaneously. Also, since I\ and Id c 
both depend linearly on KI , this calibration is indepen- 
dent of laser intensity fluctuations or reflection losses in 
optical components. For large values of ellipticity it was 
checked that this calibration method yields an agreement 
within 5% with Thorlabs PAX5710VIS, but the Thorlabs 
polarimeter cannot be reliably used to measure elliptici- 
ties smaller than 0.01 rad. 



III. JONES-MATRIX ANALYSIS 

At nearly-normal incidence, up to leading order in the 
angle of reflection 9 r , reflection coefficients for light po- 
larized parallel and perpendicular to the reflection plane 
are, respectively, 



(rg - n 2 ) 2 
(ni + n 2 ) 2 
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(6) 



This gives about 5% reflection at each optical surface for 
typical optical-window materials, in the absense of in- 
teference effects from highly parallel surfaces. Also, at 
a nearly-normal incidence, to first order in the angle of 
incidence, light is reflected without change in the polar- 
ization state. 

For the experimental setup described above (Fig. [TJ, 
in the Jones-matrix formalism, the polarization states 
for the measurement beam (ip m ) and the reference beam 
(ip r ) are given by 



ip m = M w M s M m M s M w ip , 
ijj r = M w M' m M w r)}Q, 



(7) 
(8) 



where, ipo is the Jones vector for the initial polarization 
state, and M w , M s , M m , and M' m are the Jones matri- 
ces modeling the birefringent optical-window, the sam- 
ple, the dielectric mirror, and reflection at the last op- 
tical surface of the window, respectively. In our setup 
using a dielectric mirror, M m is equal to identity as it in- 
troduces negligible ellipticity and rotation. Likewise, at 



nearly-normal incidence, M' m is equal to identity (except 
for changing overall light intensity) as well, since there 
is no measurable change in the polarization state upon 
reflection. And M s and M w are given as follows, up to a 
constant phase factor: 



M. = 
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(9) 



M w = u(e)-r Q e _^ w/2 yu- i {e\ (io) 

where (f> s and <p w are the birefringences of the sample 
and the window, respectively, 9 is the angle between the 
fast axes of the birefringent window and the sample, and 
U(9) is the rotation matrix. We assume that the sample 
birefringence has a fixed alignment, since it is controlled 
by the experimental setup. 

In the special case where the fast axis of the birefrin- 
gent window is aligned with fast axis of the sample, the 
birefringences simply add. The general case with mis- 
aligned optical axes presents a more difficult problem for 
general values of (j) a and 4> w [l3j]. However, as long as 4> s 
and 4> w are small, for tpg linearly polarized at 7r/4 rad to 
the fast axis of the sample, up to the terms linear in cf> s 
and 4> w , we obtain 

, J_ fl- j{<^ + <^(cos20 + sin20)}\ , 
V2 \l + i{4> s + (t>w(cos29-sm29)}) ' [ ' 



, 1 fl- i{^(cos 29 + sin 29)}\ 
V r ~ ^2 V + i{<Mcos26> - sin 29)} ) 



(12) 



In terms of Jones vector components, the light ellipticity 
is given by [l4[ 
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where ip x and ij) y are the complex upper and lower com- 
ponents of a Jones vector ip respectively. To first order 
in <f) s and <j) w , the ellipticities at the output are 



e m = 4>s + 4> w cos 29, 
e r = 6 W cos 29. 



(14) 
(15) 



The optical-window term appearing in the measure- 
ment beam can be independently determined from the 
light reflected from the back surface of the window. Thus, 
to first order in (j) s and <j) w , we can subtract the elliptici- 
ties of the two beams, 

e m -e r = (f) s . (16) 

In the preceding calculations, if we retain higher-order 
terms, assuming tfi s <C <ft w -C 1, the leading er- 
ror term is j4> s (f) 2 u (l — cos 49) for arbitrary orientation 
of optical-window birefringence. For example, in an ex- 
perimental setup measuring sample birefringence on the 
order of 10 -6 rad, with the optical-window birefringence 
on the order of 0.1 rad, this would lead to error of less 
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than 5 x 10 -9 rad, or a fractional error less than 0.5% of 
the sample birefringence. In contrast, for the same ex- 
perimental setup, but without the active cancellation of 
optical-window birefringence, the error would be on the 
order 0.1 rad and would overwhelm the sample signal. 

Additional considerations that may need to be taken 
into account are: 

• Rotation of polarization state of the reference beam 
introduced upon reflection at the back surface of 
the window: Since reflection coefficients as given in 
Eqs. ([S]) and © depend on the angle of reflection 
and indices of refraction, different fractions of in- 
coming light are reflected for light polarization par- 
allel and orthogonal to the fast axis of the medium, 
resulting in rotation. For an order of magnitude 
estimate of this rotation, let the fast axis of the 
medium be parallel to the plane of reflection. Let 
ni = 1 and 112 = nt for the light polarized along 
the fast axis, and ri2 — n s for the slow axis, and 
let n be average of n/ and n s . In the limit where 
(n s — rif) <Crt 2 — 1 and 6 r -C 1, when the incident 
light is polarized at 7r/4, the resulting rotation is, 

Aa=^-< (17) 
n z — 1 

The quantity (n s — rif) is related to the single-pass 
phase shift 4> w by: 

n s -n/ = (A0 U) )/(2 7 rd). (18) 

Since our modulation polarimeter setup is insen- 
sitive to rotation, we consider only the change in 
the ellipticity due to the rotation. Given a small 
rotation Aq from 7r/4, to the leading order in <f) w 
and Aa, the change in the ellipticity of the refer- 
ence beam traversing the window on the return trip 
after reflection is, 

e=Q-Aa 2 W (19) 

When the typical values (A = 632 nm, d = 1 cm, 
n = 1.5, and 9 r = 0.01 rad) are plugged in, the 
fractional change of e due to Aa is several orders 
of magnitude less than unity. 

• Ellipticity introduced into the measurement beam 
by the mirror: In the case of a dielectric mirror, 
residual stress in each dielectric layer can introduce 
additional ellipticity. At room temperature, a sep- 
arate measurement shows that the dielectric mirror 
used introduces ellipticity less than 1 mrad at inci- 
dences within 20°. 

• Spatial nonuniformity in the birefringence of win- 
dow: If the birefringence in the window is not 
uniform and the two beams pass through differ- 
ent parts of the window on the return trip, they 



may acquire different ellipticity, a systematic indis- 
tinguishable from a genuine signal. In our exper- 
iments this is the dominant source of error with 
l%-2% change in birefringence per 1 mm. 



IV. RESULTS 

The empty polarimeter has noise of 1 x 10~ 5 rad/VHz, 
but with either the window or the sample in place, the 
temperature-dependent drift in the birefringence domi- 
nates the error. Although in this experiment the ellip- 
ticity measurements of the two beams were separated in 
time by about half a minute, leading to imperfect cancel- 
lation, in principle, the drift in window birefringence can 
be canceled out with simultaneous monitoring of both 
beams. 

As a check of the experimental setup, a series of mea- 
surements was made without anything in the sample 
space. The window birefringence is varied by applying 
pressure on the fused-silica 1° optical wedge with a vise. 
In this case, we expect e s and e w , the ellipticities of mea- 
surement beam and the reference beam, respectively, to 
be equal. The obtained results (Fig. [2ji) indicate that the 
cancellation of the window birefringence is, indeed, possi- 
ble, and allow an estimate of the level to which such can- 
cellation can be done. We obtain a result consistent with 
zero difference between the two signals, and the spread of 
the points indicates that we can cancel out the window 
birefringence to about 1% of its value with the current 
setup, limited by the nonuniformity of window birefrin- 
gence and the finite beam width. 

For measurements with a sample, we take the sample 
birefringence to be the difference in ellipticity detected in 
the measurement beam when the sample is inserted and 
taken out at zero applied stress on the optical-window. 
This difference is measured at the beginning and end of 
each measurement, and the average is indicated by the 
dashed line. The two measurements also give an order- 
of-magnitude estimate of the uncertainty of the sample 
birefringence due to temperature-dependent drifts, indi- 
cated by dotted lines on the plots. For small total ellip- 
ticities, this estimate is comparable to the error of each 
data point obtained from the spread of points over a 10-s 
measurement period. 

The results for the measurement beam and the refer- 
ence beam as the stress on the optical-window is varied 
show a good agreement with the measured sample el- 
lipticity, particularly for small values of <f> s and <j) w . At 
larger values of stress on the window, the nonuniformity 
of window birefringence dominates the difference and ac- 
curate measurement of the sample becomes more diffi- 
cult. The nonuniformity of optical- window birefringence 
is measured to be about 1% change of total birefringence 
over the beam width of 1 mm for our setup. 

At larger values of 4> s and <j> w , we see the second-order 
effects due to the rotation introduced by the birefringent 
window. In principle, these effects can be compensated 
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FIG. 2: Experimental results on cancellation of window birefringence at various sample ellipticities. The sample ellipticities are as 
represented by the dashed line, with the scale on the right. The dotted lines indicate the uncertainty in determination of the sample 
birefringence-induced ellipticity change. The horizontal axis represents an increasing level of the stress on the window, and the left vertical 
axis shows the detected ellipticity. The right axis is for the difference of the ellipticity measured for the two beams. The error bars 
represent statistical errors and are obtained from the spread of the data for the measurement time of 10 s. 



for in the analysis numerically using Jones-matrix for- 
malism, provided that the azimuth a is measured in ad- 
dition to the ellipticity. Since, as shown in Eq. (p}, the 
even harmonics of u are sensitive to a, the rotation can 
be measured by using a lock-in amplifier to recover the 
second-harmonic amplitude. 

V. CONCLUSIONS 

This work demonstrates that the birefringence of a 
window can be measured using the reflection from the 
last optical surface before the sample. This measurement 



can be done simultaneously with the measurement of the 
total birefringence of the optical-window and the sam- 
ple, allowing for a real-time correction for window bire- 
fringence drifts during the experiment. This cancellation 
can effectively be done to the precision of 1% of total el- 
lipticity drift, limited by the spatial optical- window bire- 
fringence nonuniformity, for the values of optical-window 
birefringence less than 0.25 rad. 
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